This R Markdown was used to generate the figures of the publication of rMSIfragment (DOI). It serves two purposes: (1) Transparency and reproducibility of results (2) Demonstration of the main applications of rMSIfragment (for a more detailed demonstration explore demo.Rmd)
library(rMSI)
library(rMSIproc)
library(rMSIfragment)
library(ggplot2)
library(viridisLite)
rMSIfragment:::plot_figure1(res.neg.5ppm.td1,metrics[1:3],label="_NEG_D1_N")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig1_Target_Decoy_Validation_NEG_D1_N.tiff'
rMSIfragment:::plot_figure2(pks_neg,res.neg.5ppm,gt_neg,"S",steps = 1000,label="_adduct_score")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig2_Manual_Validation_adduct_score.tiff'
rMSIfragment:::plot_figure3(pks_neg,res.neg.5ppm,744.558,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
rMSIfragment:::plot_figure1(res.neg.5ppm.td2,metrics[1:3],label="_NEG_D1_N")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig1_Target_Decoy_Validation_NEG_D1_N.tiff'
rMSIfragment:::plot_figure2(pks_neg,res.neg.20ppm,gt_neg,"S",steps = 1000,label="_adduct_score")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig2_Manual_Validation_adduct_score.tiff'
rMSIfragment:::plot_figure3(pks_neg,res.neg.20ppm,744.558,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
#Compute parental
tmp<-rMSIfragment:::manual_validation(res.neg.20ppm,gt_neg)
t<-tmp$table
a<-subset(d$adducts,mode=="neg")$name
t$parental<-t$adduct%in%a
subset(t,!match&HPLC_adduct!="_")
#Quantify score for HPLC
tmp$res$match_HPLC_adduct
[1] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[25] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[49] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[73] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[97] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[121] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[145] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[169] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[193] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[217] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[241] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[265] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[289] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[313] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[337] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[361] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[385] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[409] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[433] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[457] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[481] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[505] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[529] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[553] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[577] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[601] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[625] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[649] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[673] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[697] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[721] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[745] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "[+H]+" "" "" "[+H]+" "" "[+H]+" ""
[769] "[+H]+" "" "" "" "" "" "" "" "" "" "[+H]+" "" "" "" "" "" "" "" "" "" "" "" "" ""
[793] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[817] "" "" "" "" "" "" "" "" "" "" "" "" "[+H]+" "" "" "" "" "" "" "" "[+H]+" "" "" ""
[841] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "_" "" "" "" "" "" "" "" "" ""
[865] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[889] "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "[+H]+" "" "[+H]+" "" "_" "" "" ""
[913] "" "[+H]+" "" "" "" "" "" "[+H]+" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[937] "[+H]+" "[+H]+" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" ""
[961] "" "" "_" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "_" "" "_" "" "" ""
[985] "" "" "" "" "" "" "" "" "" "[+H]+" "" "[+H]+" "" "[+H]+" "" "[+H]+"
[ reached getOption("max.print") -- omitted 8407 entries ]
i<-tmp$res$match_HPLC_adduct!=""
pval<-mean(sapply(1:100,function(j)wilcox.test(sample(tmp$res$S[!i],sum(i)), tmp$res$S[i], alternative = "two.sided")$p.value))
#roc<-td_roc(tmp$res$S,i,mz_i)
#auc<-mean(sapply(1:500,function(j)PRROC::pr.curve(scores.class1 = sample(tmp$res$S[!i],sum(i)),scores.class0=tmp$res$S[i],curve=T)$auc.integral))
pr<-PRROC::pr.curve(scores.class1 = tmp$res$S[i],scores.class0=tmp$res$S[!i],curve=T)
df<-data.frame(x=pr$curve[,1],y=pr$curve[,2])
pb<-ggplot(df,aes(x,y))+geom_line(size=2)+
xlab("Recall")+ylab("Precision")+
theme_bw()+theme(legend.justification = c(1, 0), legend.position = c(1, 0),legend.background=element_rect(size=0.2,colour=1),panel.border = element_rect(colour = "black", fill=NA), panel.grid.major = element_blank(),
panel.grid.minor = element_blank(), plot.title = element_text(hjust = 0),text=element_text(size=20)) +
ggtitle("B")#+
#scale_color_discrete(name = "Score", breaks = names,labels = paste(names," (",round(pr$auc.integral,2)," AUC)",sep = ""))
pb
pr$auc.integral
rcFinal$curve<-
rc<-PRROC::roc.curve(scores.class0 = tmp$res$S[i],scores.class1=sample(tmp$res$S[!i],sum(i)),curve=T)
df<-data.frame(x=rc$curve[,1],y=rc$curve[,2])
pc<-ggplot(df,aes(x,y))+geom_line(size=2)+
xlab("False Positive Rate")+ylab("Sensitivity")+
theme_bw()+theme(legend.justification = c(1, 0), legend.position = c(1, 0),legend.background=element_rect(size=0.2,colour=1),panel.border = element_rect(colour = "black", fill=NA), panel.grid.major = element_blank(),
panel.grid.minor = element_blank(), plot.title = element_text(hjust = 0),text=element_text(size=20)) +
ggtitle("c")#+
#scale_color_discrete(name = "Score", breaks = names,labels = paste(names," (",round(pr$auc.integral,2)," AUC)",sep = ""))
pc
rc$auc
###
rc2<-lapply(1:1000,function(j)PRROC::roc.curve(scores.class0 = tmp$res$S[i],scores.class1=sample(tmp$res$S[!i],sum(i)),curve=T))
df<-data.frame(x=unlist(lapply(rc2,function(x)x$curve[,1])),y=unlist(lapply(rc2,function(x)x$curve[,2])),id=rep(1:length(rc2),times=sapply(rc2,function(x)length(x$curve[,1]))))
ggplot(df, aes(x = x, y = y)) +
geom_line(aes(group = id), alpha = 0.5) +
geom_smooth(size = 0.8, se = F, span = 0.2)
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
#rMSIfragment:::plot_figure1(res.pos.5ppm.td1,metrics[1:3],label="_NEG_D1_N")
rMSIfragment:::plot_figure2(pks_pos,res.pos.5ppm,gt_pos,"S",steps = 1000,label="_adduct_score")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig2_Manual_Validation_adduct_score.tiff'
#rMSIfragment:::plot_figure3(pks_pos,res.pos.5ppm,744.558,3,label="_5ppm_top3")
#rMSIfragment:::plot_figure1(res.pos.5ppm.td2,metrics[1:3],label="_NEG_D1_N")
rMSIfragment:::plot_figure2(pks_pos,res.pos.20ppm,gt_pos,"S",steps = 1000,label="_adduct_score")
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig2_Manual_Validation_adduct_score.tiff'
#rMSIfragment:::plot_figure3(pks_pos,res.pos.20ppm,744.558,3,label="_5ppm_top3")
rMSIfragment:::plot_figure1(res1,metrics,label="_NEG_D1_N")
rMSIfragment:::plot_figure2(pks_neg,res,gt_neg,"adduct_score",label="_adduct_score")
rMSIfragment:::plot_figure3(pks_neg,res,744.558,3,label="_5ppm_top3")
files=list.files(path = "E:/in-source fragmentation/METASPACE/",pattern = "*.zip", full.names = TRUE,recursive = T)
pks_MS<-lapply(files,rMSIproc::LoadPeakMatrix)
csvfiles=gsub(".zip",".csv",files)
ann_MS<-lapply(csvfiles,function(x)read.csv(x,skip=2))
df<-data.frame(name=basename(files))
df$mode<-sapply(strsplit(basename(files),"_"),function(x)tolower(x[1]))
df$matrix<-sapply(strsplit(basename(files),"_"),function(x)tolower(x[2]))
df$specimen<-sapply(strsplit(basename(files),"_"),function(x)tolower(x[3]))
df$tissue<-sapply(strsplit(basename(files),"_"),function(x)tolower(x[4]))
df$analyzer<-sapply(strsplit(basename(files),"_"),function(x)tolower(x[4]))
df$mz.mean<-sapply(pks_MS,function(x)mean(x$mass))
df$mz.max<-sapply(pks_MS,function(x)max(x$mass))
df$mz.min<-sapply(pks_MS,function(x)min(x$mass))
res_MS<-lapply(seq_along(pks_MS),function(i)rMSIfragment::annotate(pks_MS[[i]],new,3,df$mode[i]))
rbind(sapply(res_MS,nrow),sapply(ann_MS,nrow),sapply(pks_MS,function(x)min(x$mass)),sapply(pks_MS,function(x)max(x$mass)))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
[1,] 237.0000 205.0000 27.0000 10021.0000 420.0000 2612.0000 596.0000 247258.0000 121.0000 2503.000 1248.0000 2203.0000 3107.0000 2406.0000 3390.0000 97.00000
[2,] 407.0000 344.0000 131.0000 319.0000 790.0000 649.0000 167.0000 889.0000 121.0000 276.000 831.0000 870.0000 2391.0000 2320.0000 3448.0000 912.00000
[3,] 202.0393 201.0242 215.0362 182.0484 299.8179 401.0702 402.7255 180.0697 206.0553 313.028 104.3449 259.1606 312.9989 365.8185 400.3417 70.06599
[4,] 887.7977 529.3906 885.8488 1903.9571 1366.2796 943.5830 1644.5788 1812.0428 783.2954 1463.211 851.8328 1281.9519 1293.2040 1293.0984 1590.7703 390.94544
i<-10
compareToMetaspace<-function(ours,metaspace,pks){
a1<-ours
a2<-subset(a1,fragmentation=="")
a3<-subset(a1,fragmentation!="")
b<-metaspace
b$mz_i<-sapply(b$mz,function(x) which.min(abs(pks$mass-x)))
matching<-paste(b$formula,b$adduct)%in%paste(a2$formula,a2$adduct)
overlappedWithFragments<-paste(b$mz_i)%in%paste(a3$mz_i)
compoundsWithFragments<-b$formula%in%a3$formula
covered<-b$formula%in%new$formula
print(nrow(ours))
return(list(covered=sum(covered)/nrow(b),matching=sum(matching)/nrow(b),overlappedWithFragments=sum(overlappedWithFragments)/nrow(b),compoundsWithFragments=sum(compoundsWithFragments)/nrow(b),v.matching=matching,v.overlappedWithFragments=overlappedWithFragments,v.compoundsWithFragments=compoundsWithFragments))
}
comp<-lapply(seq_along(pks_MS),function(i)compareToMetaspace(res_MS[[i]],ann_MS[[i]],pks_MS[[i]]))
[1] 237
[1] 205
[1] 27
[1] 10021
[1] 420
[1] 2612
[1] 596
[1] 247258
[1] 121
[1] 2503
[1] 1248
[1] 2203
[1] 3107
[1] 2406
[1] 3390
[1] 97
df$covered<-sapply(comp,function(x)x$covered)
df$matching<-sapply(comp,function(x)x$matching)
df$overlappedWithFragments<-sapply(comp,function(x)x$overlappedWithFragments)
df$compoundsWithFragments<-sapply(comp,function(x)x$compoundsWithFragments)
ggplot(df,aes(x=matching,y=covered,label=1:16))+geom_point()+geom_text(hjust=0,vjust=0)
pks_MUN<-rMSIproc::LoadPeakMatrix("E:/in-source fragmentation/MUNSTER/Munster_Au_Res_140k_10um.zip")
rMSIproc::plotPeakImage(pks_MUN,column=1000)
new<-rMSIfragment:::load_lipidMAPS("E:/in-source fragmentation/DATA/db/LMSD.sdf")
old<-rMSIfragment:::load_lipidMAPS("E:/in-source fragmentation/DATA/db/LMSDold.sdf")
tmp<-rMSIfragment:::load_decoy_non_animal(NA,"E:/in-source fragmentation/DATA/chebi/")
plot(density(tmp$exactmass))
lines(density(db_t$exactmass))
tmp2<-rMSIfragment:::adjust_db_densities(db_t,tmp)
lines(density(tmp2$exactmass))
db_d<-tmp2
save(list=c("db_t","db_d"),file="E:/in-source fragmentation/DATA/target_decoy_db_new.RData")
lapply(pks_MS,function(x)rMSIproc::plotPeakImage(x,column=900))
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
Warning in max(unlist(lapply(lstRefac[ind_offsets], function(x) { :
no non-missing arguments to max; returning -Inf
[[1]]
[[1]]$xlog
[1] FALSE
[[1]]$ylog
[1] FALSE
[[1]]$adj
[1] 0.5
[[1]]$ann
[1] TRUE
[[1]]$ask
[1] FALSE
[[1]]$bg
[1] "black"
[[1]]$bty
[1] "o"
[[1]]$cex
[1] 1
[[1]]$cex.axis
[1] 0.7
[[1]]$cex.lab
[1] 1
[[1]]$cex.main
[1] 1.2
[[1]]$cex.sub
[1] 1
[[1]]$col
[1] "white"
[[1]]$col.axis
[1] "white"
[[1]]$col.lab
[1] "white"
[[1]]$col.main
[1] "white"
[[1]]$col.sub
[1] "white"
[[1]]$crt
[1] 0
[[1]]$err
[1] 0
[[1]]$family
[1] ""
[[1]]$fg
[1] "white"
[[1]]$fig
[1] 0.000 0.875 0.000 1.000
[[1]]$fin
[1] 7.291666 4.499999
[[1]]$font
[1] 1
[[1]]$font.axis
[1] 1
[[1]]$font.lab
[1] 1
[[1]]$font.main
[1] 2
[[1]]$font.sub
[1] 1
[[1]]$lab
[1] 5 5 7
[[1]]$las
[1] 0
[[1]]$lend
[1] "round"
[[1]]$lheight
[1] 1
[[1]]$ljoin
[1] "round"
[[1]]$lmitre
[1] 10
[[1]]$lty
[1] "solid"
[[1]]$lwd
[1] 1
[[1]]$mai
[1] 0.2 0.2 0.2 0.2
[[1]]$mar
[1] 5.1 4.1 2.1 2.1
[[1]]$mex
[1] 1
[[1]]$mfcol
[1] 1 1
[[1]]$mfg
[1] 1 1 1 1
[[1]]$mfrow
[1] 1 1
[[1]]$mgp
[1] 2.0 0.5 0.5
[[1]]$mkh
[1] 0.001
[[1]]$new
[1] TRUE
[[1]]$oma
[1] 0 0 0 0
[[1]]$omd
[1] 0 1 0 1
[[1]]$omi
[1] 0 0 0 0
[[1]]$pch
[1] 1
[[1]]$pin
[1] 6.051666 3.059999
[[1]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[1]]$ps
[1] 12
[[1]]$pty
[1] "m"
[[1]]$smo
[1] 1
[[1]]$srt
[1] 0
[[1]]$tck
[1] NA
[[1]]$tcl
[1] -0.5
[[1]]$usr
[1] -32.20008 112.20008 0.00000 99.00000
[[1]]$xaxp
[1] 0 1 5
[[1]]$xaxs
[1] "r"
[[1]]$xaxt
[1] "n"
[[1]]$xpd
[1] FALSE
[[1]]$yaxp
[1] 0 1 5
[[1]]$yaxs
[1] "r"
[[1]]$yaxt
[1] "n"
[[1]]$ylbias
[1] 0.2
[[2]]
[[2]]$xlog
[1] FALSE
[[2]]$ylog
[1] FALSE
[[2]]$adj
[1] 0.5
[[2]]$ann
[1] TRUE
[[2]]$ask
[1] FALSE
[[2]]$bg
[1] "black"
[[2]]$bty
[1] "o"
[[2]]$cex
[1] 1
[[2]]$cex.axis
[1] 0.7
[[2]]$cex.lab
[1] 1
[[2]]$cex.main
[1] 1.2
[[2]]$cex.sub
[1] 1
[[2]]$col
[1] "white"
[[2]]$col.axis
[1] "white"
[[2]]$col.lab
[1] "white"
[[2]]$col.main
[1] "white"
[[2]]$col.sub
[1] "white"
[[2]]$crt
[1] 0
[[2]]$err
[1] 0
[[2]]$family
[1] ""
[[2]]$fg
[1] "white"
[[2]]$fig
[1] 0.000 0.875 0.000 1.000
[[2]]$fin
[1] 7.291666 4.499999
[[2]]$font
[1] 1
[[2]]$font.axis
[1] 1
[[2]]$font.lab
[1] 1
[[2]]$font.main
[1] 2
[[2]]$font.sub
[1] 1
[[2]]$lab
[1] 5 5 7
[[2]]$las
[1] 0
[[2]]$lend
[1] "round"
[[2]]$lheight
[1] 1
[[2]]$ljoin
[1] "round"
[[2]]$lmitre
[1] 10
[[2]]$lty
[1] "solid"
[[2]]$lwd
[1] 1
[[2]]$mai
[1] 0.2 0.2 0.2 0.2
[[2]]$mar
[1] 5.1 4.1 2.1 2.1
[[2]]$mex
[1] 1
[[2]]$mfcol
[1] 1 1
[[2]]$mfg
[1] 1 1 1 1
[[2]]$mfrow
[1] 1 1
[[2]]$mgp
[1] 2.0 0.5 0.5
[[2]]$mkh
[1] 0.001
[[2]]$new
[1] TRUE
[[2]]$oma
[1] 0 0 0 0
[[2]]$omd
[1] 0 1 0 1
[[2]]$omi
[1] 0 0 0 0
[[2]]$pch
[1] 1
[[2]]$pin
[1] 6.051666 3.059999
[[2]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[2]]$ps
[1] 12
[[2]]$pty
[1] "m"
[[2]]$smo
[1] 1
[[2]]$srt
[1] 0
[[2]]$tck
[1] NA
[[2]]$tcl
[1] -0.5
[[2]]$usr
[1] -32.20008 112.20008 0.00000 99.00000
[[2]]$xaxp
[1] 0 1 5
[[2]]$xaxs
[1] "r"
[[2]]$xaxt
[1] "n"
[[2]]$xpd
[1] FALSE
[[2]]$yaxp
[1] 0 1 5
[[2]]$yaxs
[1] "r"
[[2]]$yaxt
[1] "n"
[[2]]$ylbias
[1] 0.2
[[3]]
[[3]]$xlog
[1] FALSE
[[3]]$ylog
[1] FALSE
[[3]]$adj
[1] 0.5
[[3]]$ann
[1] TRUE
[[3]]$ask
[1] FALSE
[[3]]$bg
[1] "black"
[[3]]$bty
[1] "o"
[[3]]$cex
[1] 1
[[3]]$cex.axis
[1] 0.7
[[3]]$cex.lab
[1] 1
[[3]]$cex.main
[1] 1.2
[[3]]$cex.sub
[1] 1
[[3]]$col
[1] "white"
[[3]]$col.axis
[1] "white"
[[3]]$col.lab
[1] "white"
[[3]]$col.main
[1] "white"
[[3]]$col.sub
[1] "white"
[[3]]$crt
[1] 0
[[3]]$err
[1] 0
[[3]]$family
[1] ""
[[3]]$fg
[1] "white"
[[3]]$fig
[1] 0.000 0.875 0.000 1.000
[[3]]$fin
[1] 7.291666 4.499999
[[3]]$font
[1] 1
[[3]]$font.axis
[1] 1
[[3]]$font.lab
[1] 1
[[3]]$font.main
[1] 2
[[3]]$font.sub
[1] 1
[[3]]$lab
[1] 5 5 7
[[3]]$las
[1] 0
[[3]]$lend
[1] "round"
[[3]]$lheight
[1] 1
[[3]]$ljoin
[1] "round"
[[3]]$lmitre
[1] 10
[[3]]$lty
[1] "solid"
[[3]]$lwd
[1] 1
[[3]]$mai
[1] 0.2 0.2 0.2 0.2
[[3]]$mar
[1] 5.1 4.1 2.1 2.1
[[3]]$mex
[1] 1
[[3]]$mfcol
[1] 1 1
[[3]]$mfg
[1] 1 1 1 1
[[3]]$mfrow
[1] 1 1
[[3]]$mgp
[1] 2.0 0.5 0.5
[[3]]$mkh
[1] 0.001
[[3]]$new
[1] TRUE
[[3]]$oma
[1] 0 0 0 0
[[3]]$omd
[1] 0 1 0 1
[[3]]$omi
[1] 0 0 0 0
[[3]]$pch
[1] 1
[[3]]$pin
[1] 6.051666 3.059999
[[3]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[3]]$ps
[1] 12
[[3]]$pty
[1] "m"
[[3]]$smo
[1] 1
[[3]]$srt
[1] 0
[[3]]$tck
[1] NA
[[3]]$tcl
[1] -0.5
[[3]]$usr
[1] -29.90715 99.90715 0.00000 89.00000
[[3]]$xaxp
[1] 0 1 5
[[3]]$xaxs
[1] "r"
[[3]]$xaxt
[1] "n"
[[3]]$xpd
[1] FALSE
[[3]]$yaxp
[1] 0 1 5
[[3]]$yaxs
[1] "r"
[[3]]$yaxt
[1] "n"
[[3]]$ylbias
[1] 0.2
[[4]]
[[4]]$xlog
[1] FALSE
[[4]]$ylog
[1] FALSE
[[4]]$adj
[1] 0.5
[[4]]$ann
[1] TRUE
[[4]]$ask
[1] FALSE
[[4]]$bg
[1] "black"
[[4]]$bty
[1] "o"
[[4]]$cex
[1] 1
[[4]]$cex.axis
[1] 0.7
[[4]]$cex.lab
[1] 1
[[4]]$cex.main
[1] 1.2
[[4]]$cex.sub
[1] 1
[[4]]$col
[1] "white"
[[4]]$col.axis
[1] "white"
[[4]]$col.lab
[1] "white"
[[4]]$col.main
[1] "white"
[[4]]$col.sub
[1] "white"
[[4]]$crt
[1] 0
[[4]]$err
[1] 0
[[4]]$family
[1] ""
[[4]]$fg
[1] "white"
[[4]]$fig
[1] 0.000 0.875 0.000 1.000
[[4]]$fin
[1] 7.291666 4.499999
[[4]]$font
[1] 1
[[4]]$font.axis
[1] 1
[[4]]$font.lab
[1] 1
[[4]]$font.main
[1] 2
[[4]]$font.sub
[1] 1
[[4]]$lab
[1] 5 5 7
[[4]]$las
[1] 0
[[4]]$lend
[1] "round"
[[4]]$lheight
[1] 1
[[4]]$ljoin
[1] "round"
[[4]]$lmitre
[1] 10
[[4]]$lty
[1] "solid"
[[4]]$lwd
[1] 1
[[4]]$mai
[1] 0.2 0.2 0.2 0.2
[[4]]$mar
[1] 5.1 4.1 2.1 2.1
[[4]]$mex
[1] 1
[[4]]$mfcol
[1] 1 1
[[4]]$mfg
[1] 1 1 1 1
[[4]]$mfrow
[1] 1 1
[[4]]$mgp
[1] 2.0 0.5 0.5
[[4]]$mkh
[1] 0.001
[[4]]$new
[1] TRUE
[[4]]$oma
[1] 0 0 0 0
[[4]]$omd
[1] 0 1 0 1
[[4]]$omi
[1] 0 0 0 0
[[4]]$pch
[1] 1
[[4]]$pin
[1] 6.051666 3.059999
[[4]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[4]]$ps
[1] 12
[[4]]$pty
[1] "m"
[[4]]$smo
[1] 1
[[4]]$srt
[1] 0
[[4]]$tck
[1] NA
[[4]]$tcl
[1] -0.5
[[4]]$usr
[1] -20.05056 82.05056 0.00000 70.00000
[[4]]$xaxp
[1] 0 1 5
[[4]]$xaxs
[1] "r"
[[4]]$xaxt
[1] "n"
[[4]]$xpd
[1] FALSE
[[4]]$yaxp
[1] 0 1 5
[[4]]$yaxs
[1] "r"
[[4]]$yaxt
[1] "n"
[[4]]$ylbias
[1] 0.2
[[5]]
[[5]]$xlog
[1] FALSE
[[5]]$ylog
[1] FALSE
[[5]]$adj
[1] 0.5
[[5]]$ann
[1] TRUE
[[5]]$ask
[1] FALSE
[[5]]$bg
[1] "black"
[[5]]$bty
[1] "o"
[[5]]$cex
[1] 1
[[5]]$cex.axis
[1] 0.7
[[5]]$cex.lab
[1] 1
[[5]]$cex.main
[1] 1.2
[[5]]$cex.sub
[1] 1
[[5]]$col
[1] "white"
[[5]]$col.axis
[1] "white"
[[5]]$col.lab
[1] "white"
[[5]]$col.main
[1] "white"
[[5]]$col.sub
[1] "white"
[[5]]$crt
[1] 0
[[5]]$err
[1] 0
[[5]]$family
[1] ""
[[5]]$fg
[1] "white"
[[5]]$fig
[1] 0.000 0.875 0.000 1.000
[[5]]$fin
[1] 7.291666 4.499999
[[5]]$font
[1] 1
[[5]]$font.axis
[1] 1
[[5]]$font.lab
[1] 1
[[5]]$font.main
[1] 2
[[5]]$font.sub
[1] 1
[[5]]$lab
[1] 5 5 7
[[5]]$las
[1] 0
[[5]]$lend
[1] "round"
[[5]]$lheight
[1] 1
[[5]]$ljoin
[1] "round"
[[5]]$lmitre
[1] 10
[[5]]$lty
[1] "solid"
[[5]]$lwd
[1] 1
[[5]]$mai
[1] 0.2 0.2 0.2 0.2
[[5]]$mar
[1] 5.1 4.1 2.1 2.1
[[5]]$mex
[1] 1
[[5]]$mfcol
[1] 1 1
[[5]]$mfg
[1] 1 1 1 1
[[5]]$mfrow
[1] 1 1
[[5]]$mgp
[1] 2.0 0.5 0.5
[[5]]$mkh
[1] 0.001
[[5]]$new
[1] TRUE
[[5]]$oma
[1] 0 0 0 0
[[5]]$omd
[1] 0 1 0 1
[[5]]$omi
[1] 0 0 0 0
[[5]]$pch
[1] 1
[[5]]$pin
[1] 6.051666 3.059999
[[5]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[5]]$ps
[1] 12
[[5]]$pty
[1] "m"
[[5]]$smo
[1] 1
[[5]]$srt
[1] 0
[[5]]$tck
[1] NA
[[5]]$tcl
[1] -0.5
[[5]]$usr
[1] -29.90715 99.90715 0.00000 89.00000
[[5]]$xaxp
[1] 0 1 5
[[5]]$xaxs
[1] "r"
[[5]]$xaxt
[1] "n"
[[5]]$xpd
[1] FALSE
[[5]]$yaxp
[1] 0 1 5
[[5]]$yaxs
[1] "r"
[[5]]$yaxt
[1] "n"
[[5]]$ylbias
[1] 0.2
[[6]]
[[6]]$xlog
[1] FALSE
[[6]]$ylog
[1] FALSE
[[6]]$adj
[1] 0.5
[[6]]$ann
[1] TRUE
[[6]]$ask
[1] FALSE
[[6]]$bg
[1] "black"
[[6]]$bty
[1] "o"
[[6]]$cex
[1] 1
[[6]]$cex.axis
[1] 0.7
[[6]]$cex.lab
[1] 1
[[6]]$cex.main
[1] 1.2
[[6]]$cex.sub
[1] 1
[[6]]$col
[1] "white"
[[6]]$col.axis
[1] "white"
[[6]]$col.lab
[1] "white"
[[6]]$col.main
[1] "white"
[[6]]$col.sub
[1] "white"
[[6]]$crt
[1] 0
[[6]]$err
[1] 0
[[6]]$family
[1] ""
[[6]]$fg
[1] "white"
[[6]]$fig
[1] 0.000 0.875 0.000 1.000
[[6]]$fin
[1] 7.291666 4.499999
[[6]]$font
[1] 1
[[6]]$font.axis
[1] 1
[[6]]$font.lab
[1] 1
[[6]]$font.main
[1] 2
[[6]]$font.sub
[1] 1
[[6]]$lab
[1] 5 5 7
[[6]]$las
[1] 0
[[6]]$lend
[1] "round"
[[6]]$lheight
[1] 1
[[6]]$ljoin
[1] "round"
[[6]]$lmitre
[1] 10
[[6]]$lty
[1] "solid"
[[6]]$lwd
[1] 1
[[6]]$mai
[1] 0.2 0.2 0.2 0.2
[[6]]$mar
[1] 5.1 4.1 2.1 2.1
[[6]]$mex
[1] 1
[[6]]$mfcol
[1] 1 1
[[6]]$mfg
[1] 1 1 1 1
[[6]]$mfrow
[1] 1 1
[[6]]$mgp
[1] 2.0 0.5 0.5
[[6]]$mkh
[1] 0.001
[[6]]$new
[1] TRUE
[[6]]$oma
[1] 0 0 0 0
[[6]]$omd
[1] 0 1 0 1
[[6]]$omi
[1] 0 0 0 0
[[6]]$pch
[1] 1
[[6]]$pin
[1] 6.051666 3.059999
[[6]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[6]]$ps
[1] 12
[[6]]$pty
[1] "m"
[[6]]$smo
[1] 1
[[6]]$srt
[1] 0
[[6]]$tck
[1] NA
[[6]]$tcl
[1] -0.5
[[6]]$usr
[1] 0.000000 167.000000 -1.247164 113.247164
[[6]]$xaxp
[1] 0 1 5
[[6]]$xaxs
[1] "r"
[[6]]$xaxt
[1] "n"
[[6]]$xpd
[1] FALSE
[[6]]$yaxp
[1] 0 1 5
[[6]]$yaxs
[1] "r"
[[6]]$yaxt
[1] "n"
[[6]]$ylbias
[1] 0.2
[[7]]
[[7]]$xlog
[1] FALSE
[[7]]$ylog
[1] FALSE
[[7]]$adj
[1] 0.5
[[7]]$ann
[1] TRUE
[[7]]$ask
[1] FALSE
[[7]]$bg
[1] "black"
[[7]]$bty
[1] "o"
[[7]]$cex
[1] 1
[[7]]$cex.axis
[1] 0.7
[[7]]$cex.lab
[1] 1
[[7]]$cex.main
[1] 1.2
[[7]]$cex.sub
[1] 1
[[7]]$col
[1] "white"
[[7]]$col.axis
[1] "white"
[[7]]$col.lab
[1] "white"
[[7]]$col.main
[1] "white"
[[7]]$col.sub
[1] "white"
[[7]]$crt
[1] 0
[[7]]$err
[1] 0
[[7]]$family
[1] ""
[[7]]$fg
[1] "white"
[[7]]$fig
[1] 0.000 0.875 0.000 1.000
[[7]]$fin
[1] 7.291666 4.499999
[[7]]$font
[1] 1
[[7]]$font.axis
[1] 1
[[7]]$font.lab
[1] 1
[[7]]$font.main
[1] 2
[[7]]$font.sub
[1] 1
[[7]]$lab
[1] 5 5 7
[[7]]$las
[1] 0
[[7]]$lend
[1] "round"
[[7]]$lheight
[1] 1
[[7]]$ljoin
[1] "round"
[[7]]$lmitre
[1] 10
[[7]]$lty
[1] "solid"
[[7]]$lwd
[1] 1
[[7]]$mai
[1] 0.2 0.2 0.2 0.2
[[7]]$mar
[1] 5.1 4.1 2.1 2.1
[[7]]$mex
[1] 1
[[7]]$mfcol
[1] 1 1
[[7]]$mfg
[1] 1 1 1 1
[[7]]$mfrow
[1] 1 1
[[7]]$mgp
[1] 2.0 0.5 0.5
[[7]]$mkh
[1] 0.001
[[7]]$new
[1] TRUE
[[7]]$oma
[1] 0 0 0 0
[[7]]$omd
[1] 0 1 0 1
[[7]]$omi
[1] 0 0 0 0
[[7]]$pch
[1] 1
[[7]]$pin
[1] 6.051666 3.059999
[[7]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[7]]$ps
[1] 12
[[7]]$pty
[1] "m"
[[7]]$smo
[1] 1
[[7]]$srt
[1] 0
[[7]]$tck
[1] NA
[[7]]$tcl
[1] -0.5
[[7]]$usr
[1] -4.878596e+01 1.247860e+02 7.105427e-15 1.190000e+02
[[7]]$xaxp
[1] 0 1 5
[[7]]$xaxs
[1] "r"
[[7]]$xaxt
[1] "n"
[[7]]$xpd
[1] FALSE
[[7]]$yaxp
[1] 0 1 5
[[7]]$yaxs
[1] "r"
[[7]]$yaxt
[1] "n"
[[7]]$ylbias
[1] 0.2
[[8]]
[[8]]$xlog
[1] FALSE
[[8]]$ylog
[1] FALSE
[[8]]$adj
[1] 0.5
[[8]]$ann
[1] TRUE
[[8]]$ask
[1] FALSE
[[8]]$bg
[1] "black"
[[8]]$bty
[1] "o"
[[8]]$cex
[1] 1
[[8]]$cex.axis
[1] 0.7
[[8]]$cex.lab
[1] 1
[[8]]$cex.main
[1] 1.2
[[8]]$cex.sub
[1] 1
[[8]]$col
[1] "white"
[[8]]$col.axis
[1] "white"
[[8]]$col.lab
[1] "white"
[[8]]$col.main
[1] "white"
[[8]]$col.sub
[1] "white"
[[8]]$crt
[1] 0
[[8]]$err
[1] 0
[[8]]$family
[1] ""
[[8]]$fg
[1] "white"
[[8]]$fig
[1] 0.000 0.875 0.000 1.000
[[8]]$fin
[1] 7.291666 4.499999
[[8]]$font
[1] 1
[[8]]$font.axis
[1] 1
[[8]]$font.lab
[1] 1
[[8]]$font.main
[1] 2
[[8]]$font.sub
[1] 1
[[8]]$lab
[1] 5 5 7
[[8]]$las
[1] 0
[[8]]$lend
[1] "round"
[[8]]$lheight
[1] 1
[[8]]$ljoin
[1] "round"
[[8]]$lmitre
[1] 10
[[8]]$lty
[1] "solid"
[[8]]$lwd
[1] 1
[[8]]$mai
[1] 0.2 0.2 0.2 0.2
[[8]]$mar
[1] 5.1 4.1 2.1 2.1
[[8]]$mex
[1] 1
[[8]]$mfcol
[1] 1 1
[[8]]$mfg
[1] 1 1 1 1
[[8]]$mfrow
[1] 1 1
[[8]]$mgp
[1] 2.0 0.5 0.5
[[8]]$mkh
[1] 0.001
[[8]]$new
[1] TRUE
[[8]]$oma
[1] 0 0 0 0
[[8]]$omd
[1] 0 1 0 1
[[8]]$omi
[1] 0 0 0 0
[[8]]$pch
[1] 1
[[8]]$pin
[1] 6.051666 3.059999
[[8]]$plt
[1] 0.08502858 0.91497142 0.16000004 0.83999996
[[8]]$ps
[1] 12
[[8]]$pty
[1] "m"
[[8]]$smo
[1] 1
[[8]]$srt
[1] 0
[[8]]$tck
[1] NA
[[8]]$tcl
[1] -0.5
[[8]]$usr
[1] -5.424455e+01 1.222445e+02 7.105427e-15 1.210000e+02
[[8]]$xaxp
[1] 0 1 5
[[8]]$xaxs
[1] "r"
[[8]]$xaxt
[1] "n"
[[8]]$xpd
[1] FALSE
[[8]]$yaxp
[1] 0 1 5
[[8]]$yaxs
[1] "r"
[[8]]$yaxt
[1] "n"
[[8]]$ylbias
[1] 0.2
1
[1] 1
matched_us<-matched_us[order(matched_us$S),]
matched_us<-matched_us[order(matched_us$S),]
rMSIfragment:::plot_figure3(pks_MS[[4]],res_MS[[4]],772.5864,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
rMSIfragment:::plot_figure3(pks_MS[[4]],res_MS[[4]],744.5553,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
rMSIfragment:::plot_figure3(pks_MS[[4]],res_MS[[4]],887.5661,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
rMSIfragment:::plot_figure3(pks_MS[[4]],res_MS[[4]],788.5452,3,label="_5ppm_top3")
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
Warning in grDevices::dev.off() :
unable to open TIFF file '/home/gbaquer/msidata/1. In-source Fragmentation/1.6. Paper/Figures//Fig3_Example_Annotations_5ppm_top3.tiff'
them[which.min(abs(them$mz-887.57)),]
them[which.min(abs(them$mz-836.53)),]
them[which.min(abs(them$mz-749.51)),]
files_imzML=list.files(path = "E:/in-source fragmentation/METASPACE/new3/",pattern = "*.imzML", full.names = TRUE,recursive = T)
files_zip=gsub("*.imzML",".zip",files_imzML)
for(i in seq_along(files_imzML)){
tmp<-rMSIworkflows:::loadBinnedCentroidData(files_imzML[i],3)
rMSIproc::StorePeakMatrix(files_zip[i],tmp)
}
Checking binary data checksum using MD5 key... OK
Reading the .ibd file(binary)...
|
| | 0%
|
|= | 0%
|
|= | 1%
|
|== | 1%
|
|=== | 1%
|
|=== | 2%
|
|==== | 2%
|
|===== | 2%
|
|===== | 3%
|
|====== | 3%
|
|======= | 3%
|
|======= | 4%
|
|======== | 4%
|
|========= | 4%
|
|========= | 5%
|
|========== | 5%
|
|=========== | 6%
|
|============ | 6%
|
|============ | 7%
|
|============= | 7%
|
|============== | 7%
|
|============== | 8%
|
|=============== | 8%
|
|================ | 8%
|
|================ | 9%
|
|================= | 9%
|
|================== | 9%
|
|================== | 10%
|
|=================== | 10%
|
|==================== | 10%
|
|==================== | 11%
|
|===================== | 11%
|
|====================== | 11%
|
|====================== | 12%
|
|======================= | 12%
|
|======================== | 12%
|
|======================== | 13%
|
|========================= | 13%
|
|========================== | 13%
|
|========================== | 14%
|
|=========================== | 14%
|
|============================ | 14%
|
|============================ | 15%
|
|============================= | 15%
|
|============================== | 15%
|
|============================== | 16%
|
|=============================== | 16%
|
|================================ | 17%
|
|================================= | 17%
|
|================================= | 18%
|
|================================== | 18%
|
|=================================== | 18%
|
|=================================== | 19%
|
|==================================== | 19%
|
|===================================== | 19%
|
|===================================== | 20%
|
|====================================== | 20%
|
|======================================= | 20%
|
|======================================= | 21%
|
|======================================== | 21%
|
|========================================= | 21%
|
|========================================= | 22%
|
|========================================== | 22%
|
|=========================================== | 22%
|
|=========================================== | 23%
|
|============================================ | 23%
|
|============================================= | 23%
|
|============================================= | 24%
|
|============================================== | 24%
|
|=============================================== | 24%
|
|=============================================== | 25%
|
|================================================ | 25%
|
|================================================= | 25%
|
|================================================= | 26%
|
|================================================== | 26%
|
|=================================================== | 26%
|
|=================================================== | 27%
|
|==================================================== | 27%
|
|===================================================== | 28%
|
|====================================================== | 28%
|
|====================================================== | 29%
|
|======================================================= | 29%
|
|======================================================== | 29%
|
|======================================================== | 30%
|
|========================================================= | 30%
|
|========================================================== | 30%
|
|========================================================== | 31%
|
|=========================================================== | 31%
|
|============================================================ | 31%
|
|============================================================ | 32%
|
|============================================================= | 32%
|
|============================================================== | 32%
|
|============================================================== | 33%
|
|=============================================================== | 33%
|
|================================================================ | 33%
|
|================================================================ | 34%
|
|================================================================= | 34%
|
|================================================================== | 34%
|
|================================================================== | 35%
|
|=================================================================== | 35%
|
|==================================================================== | 35%
|
|==================================================================== | 36%
|
|===================================================================== | 36%
|
|====================================================================== | 36%
|
|====================================================================== | 37%
|
|======================================================================= | 37%
|
|======================================================================== | 37%
|
|======================================================================== | 38%
|
|========================================================================= | 38%
|
|========================================================================== | 38%
|
|========================================================================== | 39%
|
|=========================================================================== | 39%
|
|=========================================================================== | 40%
|
|============================================================================ | 40%
|
|============================================================================= | 40%
|
|============================================================================= | 41%
|
|============================================================================== | 41%
|
|=============================================================================== | 41%
|
|=============================================================================== | 42%
|
|================================================================================ | 42%
|
|================================================================================= | 42%
|
|================================================================================= | 43%
|
|================================================================================== | 43%
|
|=================================================================================== | 43%
|
|=================================================================================== | 44%
|
|==================================================================================== | 44%
|
|===================================================================================== | 44%
|
|===================================================================================== | 45%
|
|====================================================================================== | 45%
|
|======================================================================================= | 45%
|
|======================================================================================= | 46%
|
|======================================================================================== | 46%
|
|========================================================================================= | 46%
|
|========================================================================================= | 47%
|
|========================================================================================== | 47%
|
|=========================================================================================== | 47%
|
|=========================================================================================== | 48%
|
|============================================================================================ | 48%
|
|============================================================================================= | 48%
|
|============================================================================================= | 49%
|
|============================================================================================== | 49%
|
|=============================================================================================== | 49%
|
|=============================================================================================== | 50%
|
|================================================================================================ | 50%
|
|================================================================================================ | 51%
|
|================================================================================================= | 51%
|
|================================================================================================== | 51%
|
|================================================================================================== | 52%
|
|=================================================================================================== | 52%
|
|==================================================================================================== | 52%
|
|==================================================================================================== | 53%
|
|===================================================================================================== | 53%
|
|====================================================================================================== | 53%
|
|====================================================================================================== | 54%
|
|======================================================================================================= | 54%
|
|======================================================================================================== | 54%
|
|======================================================================================================== | 55%
|
|========================================================================================================= | 55%
|
|========================================================================================================== | 55%
|
|========================================================================================================== | 56%
|
|=========================================================================================================== | 56%
|
|============================================================================================================ | 56%
|
|============================================================================================================ | 57%
|
|============================================================================================================= | 57%
|
|============================================================================================================== | 57%
|
|============================================================================================================== | 58%
|
|=============================================================================================================== | 58%
|
|================================================================================================================ | 58%
|
|================================================================================================================ | 59%
|
|================================================================================================================= | 59%
|
|================================================================================================================== | 59%
|
|================================================================================================================== | 60%
|
|=================================================================================================================== | 60%
|
|==================================================================================================================== | 60%
|
|==================================================================================================================== | 61%
|
|===================================================================================================================== | 61%
|
|===================================================================================================================== | 62%
|
|====================================================================================================================== | 62%
|
|======================================================================================================================= | 62%
|
|======================================================================================================================= | 63%
|
|======================================================================================================================== | 63%
|
|========================================================================================================================= | 63%
|
|========================================================================================================================= | 64%
|
|========================================================================================================================== | 64%
|
|=========================================================================================================================== | 64%
|
|=========================================================================================================================== | 65%
|
|============================================================================================================================ | 65%
|
|============================================================================================================================= | 65%
|
|============================================================================================================================= | 66%
|
|============================================================================================================================== | 66%
|
|=============================================================================================================================== | 66%
|
|=============================================================================================================================== | 67%
|
|================================================================================================================================ | 67%
|
|================================================================================================================================= | 67%
|
|================================================================================================================================= | 68%
|
|================================================================================================================================== | 68%
|
|=================================================================================================================================== | 68%
|
|=================================================================================================================================== | 69%
|
|==================================================================================================================================== | 69%
|
|===================================================================================================================================== | 69%
|
|===================================================================================================================================== | 70%
|
|====================================================================================================================================== | 70%
|
|======================================================================================================================================= | 70%
|
|======================================================================================================================================= | 71%
|
|======================================================================================================================================== | 71%
|
|========================================================================================================================================= | 71%
|
|========================================================================================================================================= | 72%
|
|========================================================================================================================================== | 72%
|
|=========================================================================================================================================== | 73%
|
|============================================================================================================================================ | 73%
|
|============================================================================================================================================ | 74%
|
|============================================================================================================================================= | 74%
|
|============================================================================================================================================== | 74%
|
|============================================================================================================================================== | 75%
|
|=============================================================================================================================================== | 75%
|
|================================================================================================================================================ | 75%
|
|================================================================================================================================================ | 76%
|
|================================================================================================================================================= | 76%
|
|================================================================================================================================================== | 76%
|
|================================================================================================================================================== | 77%
|
|=================================================================================================================================================== | 77%
|
|==================================================================================================================================================== | 77%
|
|==================================================================================================================================================== | 78%
|
|===================================================================================================================================================== | 78%
|
|====================================================================================================================================================== | 78%
|
|====================================================================================================================================================== | 79%
|
|======================================================================================================================================================= | 79%
|
|======================================================================================================================================================== | 79%
|
|======================================================================================================================================================== | 80%
|
|========================================================================================================================================================= | 80%
|
|========================================================================================================================================================== | 80%
|
|========================================================================================================================================================== | 81%
|
|=========================================================================================================================================================== | 81%
|
|============================================================================================================================================================ | 81%
|
|============================================================================================================================================================ | 82%
|
|============================================================================================================================================================= | 82%
|
|============================================================================================================================================================== | 82%
|
|============================================================================================================================================================== | 83%
|
|=============================================================================================================================================================== | 83%
|
|================================================================================================================================================================ | 84%
|
|================================================================================================================================================================= | 84%
|